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Abstract 

We introduce quantum fluctuations into the simulated annealing process 
of optimization problems, aiming at faster convergence to the optimal state. 
Quantum fluctuations cause transitions between states and thus play the same 
role as thermal fluctuations in the conventional approach. The idea is tested 
by the transverse Ising model, in which the transverse field is a function of 
time similar to the temperature in the conventional method. The goal is to 
find the ground state of the diagonal part of the Hamiltonian with high accu- 
racy as quickly as possible. We have solved the time-dependent Schrodinger 
equation numerically for small size systems with various exchange interac- 
tions. Comparison with the results of the corresponding classical (thermal) 
method reveals that the quantum annealing leads to the ground state with 
much larger probability in almost all cases if we use the same annealing sched- 
ule. 
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I. INTRODUCTION 



The technique of simulated annealing (SA) was first proposed by Kirkpatrick et al. (T| as 
a general method to solve optimization problems. The idea is to use thermal fluctuations to 
allow the system to escape from local minima of the cost function so that the system reaches 
the global minimum under an appropriate annealing schedule (the rate of decrease of tem- 
perature). If the temperature is decreased too quickly, the system may become trapped in 
a local minimum. Too slow annealing, on the other hand, is practically useless although 
such a process would certainly bring the system to the global minimum. Geman and Ge- 
man proved a theorem on the annealing schedule for a generic problem of combinatorial 
optimization ||. They showed that any system reaches the global minimum of the cost 
function asymptotically if the temperature is decreased as T = c/ log t or slower, where c is 
a constant determined by the system size and other structures of the cost function. This 
bound on the annealing schedule may be the optimal one under generic conditions although 
faster decrease of the temperature often gives satisfactory results in practical applications 
for many systems. 

Thermal fluctuations were introduced in the optimization problem so that transitions 
between states take place in the process of search for the global minimum among many states. 
Thus there seems to be no reasons to avoid use of other mechanisms for state transitions 
if these mechanisms may lead to better convergence properties. One such possibility is the 
generalized transition probability of Tsallis [[J , which is a generalization of the conventional 
Boltzmann-type transition probability appearing in the master equation and thus used in 
Monte Carlo simulations (see also [|J). In the present paper we seek another possibility of 
making use of quantum tunneling processes for state transitions, which we call quantum 
annealing (QA). In particular we would like to learn how effectively quantum tunneling 
processes possibly lead to the global minimum in comparison to temperature-driven processes 
used in the conventional method of SA. 

In the virtual absence of previous studies along such a line of consideration, it seems 
better to focus our attention on a specific model system, rather than to develop a general 
argument, to gain an insight into the role of quantum fluctuations in the situation of opti- 
mization problem. Quantum effects have been found to play a very similar role to thermal 
fluctuations in the Hopfield model in a transverse field in thermal equilibrium ||. This 
observation motivates us to investigate dynamical properties of the Ising model under quan- 
tum fluctuations in the form of a transverse field. We therefore discuss in this paper the 
transverse Ising model with a variety of exchange interactions. The transverse field controls 
the rate of transition between states and thus plays the same role as the temperature does in 
SA. We assume that the system has no thermal fluctuations in the QA context and the term 
"ground state" refers to the lowest-energy state of the Hamiltonian without the transverse 
field term. 

Static properties of the transverse Ising model have been investigated quite extensively 
for many years ||. There have however been very few studies on the dynamical behavior of 
the Ising model with a transverse field. We refer to the work by Sato et al. who carried out 
quantum Monte Carlo simulations of the two dimensional gaussian spin glass model in an 
infinitesimal transverse field, showing a reasonably fast approach to the ground state J7|. 

We present here a new point of view by comparing the efficiency of QA directly with 
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that of classical SA in reaching the ground state. We solve the time-dependent Schrodinger 
equation and the classical master equation numerically for small-size systems with the same 
exchange interactions under the same annealing schedules. Calculations of probabilities that 
the system is in the ground state at each time for both classical and quantum cases give 
important implications on the relative efficiency of the two approaches. 

In the next section we explain the model and define the measure of closeness of the 
system in QA to the desired ground state. This measure is compared with the corresponding 
classical probability that the system is in the ground state, the definition of which is also 
given there. In Sec. [Ill] numerical results for QA and SA are shown for various types of 
annealing schedules and interactions. The data suggest that QA generally gives a larger 
probability to lead to the ground state than SA under the same conditions on the annealing 
schedule and interactions. Section [TV] deals with the analytical solutions for the one-spin 
case which turns out to be quite non-trivial. Explicit solutions yield very useful information 
to clarify several subtle aspects of the problem. The final section is devoted to summary 
and discussions. 



II. TRANSVERSE ISING MODEL 

Let us consider the following Ising model with longitudinal and transverse fields: 

m = -£ j^to] -h^i-mT. < (!) 

ij i i 

= H -r(t)J2af, (2) 



where the type of interactions will be specified later. The term of longitudinal field was 
introduced to remove the trivial degeneracy in the exchange interaction term coming from 
the overall up-down symmetry that effectively reduces the available phase space by half. 
The r(i)-term causes quantum tunneling between various classical states (the eigenstates of 
the classical part TCo)- By decreasing the amplitude T(t) of the transverse field from a very 
large value to zero, we hopefully drive the system into the optimal state, the ground state 
ofH . 

The natural dynamics of the present system is provided by the Schrodinger equation 

= (3) 

We solve this time-dependent Schrodinger equation numerically for small-size systems. The 
representation to diagonalize Tio (the ^-representation) will be used throughout the paper. 
The corresponding classical SA process is described by the master equation 

^ = E/vm (4) 

where P%{t) represents the probability that the system is in the zth state. We consider 
single-spin flip processes with the transition matrix elements given as 



3 




{1 + exp((Ei — Ej)/T(t))} 1 (single-spin difference) 
(otherwise) 



(5) 



In SA, the temperature T(t) is first set to a very large value and then is gradually 
decreased to zero. The corresponding process in QA should be to change T(t) from a very 
large value to zero. The reason is that the high-temperature state in SA is a mixture of 
all possible states with almost equal probabilities, and the corresponding state in QA is the 
linear combination of all states with equal amplitude in the ^-representation, which is the 
lowest eigenstate of the Hamiltonian (|I|) for very large T. The low-temperature state after a 
successful SA is the ground state of 7i , which should also be the eigenstate of H(t) as T(t) 
is reduced to zero sufficiently slowly in QA. Another justification of identification of T and 
T comes from the fact that the T = phase diagram of the Hopfield model in a transverse 
field has almost the same structure as the equilibrium phase diagram of the conventional 
Hopfield model at finite temperature if we identify the temperature axis of the latter phase 
diagram with the T axis in the former |J. We therefore change T(t) in QA and T(t) in SA 
from infinity to zero with the same functional forms T(t) = T(t) = c/t, c/\/i, cj log(t + 1) 
(t : — > oo) or —ct (t : — oo — > 0). The reason for choosing these functional forms are that 
either they allow for analytical solutions in the single-spin case as shown in Sec. [TV] or for 
comparison with the Geman-Geman bound mentioned in Sec. |I|. 

To compare the performance of the two methods QA and SA, we calculate the probabili- 
ties -Pqa(^) = | (gl^it)) | 2 for QA and Psa(^) = Pg{t) for SA, where P g (t) is the probability to 
find the system in the ground state at time t in SA and \g) is the ground-state wave function 
of 7^o- Note that we treat only small-size systems (the number of spins N = 8) and thus 
the ground state can be picked out explicitly. In the ideal situation Pqa(£) and Psa(^) will 
be very small initially and increase towards 1 as t —>■ oo. 

It is useful to introduce another set of quantities P| A (T) and Pq A (T). The former 
is the Boltzmann factor of the ground state of 7io at temperature T while the latter is 
defined as | (g\ipr) f, where the wave function ipr is the lowest-energy stationary state of 
the full Hamiltonian (|l|) for a given fixed value of T. In the quasi-static limit, the system 
follows equilibrium in SA and thus -Psa(^) = PsA.(T(t)). Correspondingly for QA, -Pqa(^) = 
Pq A (r(t)) when T(t) changes sufficiently slowly. Thus the differences between both sides of 
these two equations give measures how closely the system follows quasi-static states during 
dynamical process of annealing. 



We now present numerical results on Psa and Pqa for various types of exchange inter- 
actions and transverse fields. All calculations were performed with a constant longitudinal 
field h = 0.1 to remove trivial degeneracy. 



Let us first discuss the ferromagnetic Ising model with J = const for all pairs of spins. 
Figure [l] shows the overlaps for the case of T(t) = T(t) = 3/log(t + 1). It is seen that 



III. NUMERICAL RESULTS 



A. Ferromagnetic Model 
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both QA and SA follow stationary (equilibrium) states during dynamical processes rather 
accurately. In SA the theorem of Geman and Geman [|2j guarantees that the annealing 
schedule T(t) = c/log(l + t) assures convergence to the ground state (Psa — > 1 in our 
notation) if c is adjusted appropriately. Our choice c = 3 is somewhat arbitrary but the 
tendency is clear for Psa — ► 1 as t — > oo, which is also clear from approximate satisfaction of 
the quasi-equilibrium condition Psa(^) = -P|a(^(^))- Although there are no mathematically 
rigorous arguments for QA corresponding to the Geman-Geman bound, the numerical data 
indicate convergence to the ground state under the annealing schedule T(t) = 3/log(t + 1) 
at least for the ferromagnetic system. It should be remembered that the unit of time is 
arbitrary since we have set h — 1 in the Schrddinger equation (|3]) and the unit of time r = 1 
in the master equation (|j). Thus the fact that the curves for QA in Fig. [l| lie below those 
for SA at any given time does not have any positive significance. 

If we decrease the transverse field and the temperature faster, T(t) = T(t) = 3/ \ft, there 
appears a qualitative difference between QA and SA as shown in Fig. ^. The quantum 
method clearly gives better convergence to the ground state while the classical counterpart 
gets stuck in a local minimum with a non-negligible probability. To see the rate of approach 
of Pqa to 1, we have plotted 1 — Pqa in a log- log scale in Fig. |^. It is seen that 1 — Pqa 
behaves as const/t in the time region between 100 and 1000. 

By a still faster annealing schedule T{t) = T(t) = 3/i, the system becomes trapped in 
intermediate states both in QA and SA as seen in Fig. f|. 

B. Frustrated Model 

We next analyze the interesting case of a frustrated system shown in Fig. |5[ The full 
lines indicate ferromagnetic interactions while the broken line is for an antiferromagnetic 
interaction with the same absolute value as the ferromagnetic ones. If the temperature is 
very high in the classical case, the spins 4 and 5 are changing their states very rapidly and 
hence the effective interaction between spins 3 and 6 via spins 4 and 5 will be negligibly 
small. Thus the direct antiferromagnetic interactions between spins 3 and 6 is expected to 
dominate the correlation of these spins, which is clearly observed in Fig. |6] as the negative 
value of the thermodynamic correlation function (cr|<7g) c in the high-temperature side. At 
low temperatures, on the other hand, the spins 4 and 5 tend to be fixed in some definite 
direction and consequently the effective ferromagnetic interactions between spins 3 and 6 are 
roughly twice as large as the direct antiferromagnetic interaction. This argument is justified 
by the positive value of the correlation function at low temperatures in Fig. || Therefore the 
spins 3 and 6 must change their relative orientation at some intermediate temperature. This 
means that the free-energy landscape goes under significant restructuring as the temperature 
is decreased and therefore the annealing process should be performed with sufficient care. 

If the transverse field in QA plays a similar role to the temperature in SA, we expect 
similar dependence of the correlation function (cr|o"|) g on the transverse field T. Here the 
expectation value is evaluated by the stationary eigenfunction of the full Hamiltonian ([[]) 
with the lowest eigenvalue at a given T. The broken curve in Fig. ^ clearly supports this 
idea. We therefore expect that the frustrated system of Fig. ^] is a good test ground for 
comparison of QA and SA in the situation with a significant change of spin configurations 
in the dynamical process of annealing. 
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The results are shown in Fig. |7| for the annealing schedule T(t) = T(t) = 3/y/i. The 
time scale is normalized so that both classical and quantum correlation functions vanish at 
t = 1 (see Fig. H). The tendency is clear that QA is better suited for ground-state search in 
the present system. 



C. Random Interaction Model 



The third and final example is the Sherrington-Kirkpatrick (SK) model of spin glasses 
||. Interactions exist between all pairs of spins and are chosen from a Gaussian distribution 
with vanishing mean and variance 1/N (N = 8 in our case). Figure |8| shows a typical result 
on the time evolution of the probabilities under the annealing schedule T(t) = T(t) = 3/ 
We have checked several realizations of exchange interactions under the same distribution 
function and have found that the results are qualitatively the same. Figure || again suggests 
that QA is better suited than SA for the present optimization problem. 



IV. SOLUTION OF THE SINGLE-SPIN PROBLEM 

It is possible to solve the time-dependent Schrodinger equation explicitly when the prob- 
lem involves only a single spin and the functional form of the transverse field is T(t) = —ct, cj t 
or c/\/i. We note that the single-spin problem is trivial in SA because there are only two 
states involved (up and down) and thus there are no local minima. This does not mean that 
the same single-spin problem is also trivial in the quantum mechanical version. In QA with 
a single spin, the transition between the two states is caused by a finite transverse field. 
The system goes through tunneling processes to reach the other state, and an approximate 
annealing schedule is essential to reach the ground state. On the other hand, in SA, the 
transition from the higher state to the lower state takes place even at T = and thus the 
system always reaches the ground state. 

Let us first discuss the case of T(t) = —ct with t changing from — oo to 0. This is the 
well-known Landau-Zener model and the explicit solution of the time-dependent Schrodinger 
equation is available in the literature [p|-[l3||. With the notation a(t) = (+\tp(t)) and b(t) = 
(—\ijj{t)) and the initial condition o(— oo) = b(— oo) = l/y/2 (the lowest eigenstate), the 
solution for b(t) is found to be (see Appendix) 

he~* h2 / Sc f 2ct + h . 
b(t) = j^D^i-tz) 

-'-r^ ^P \ 2 ( /.:)). (6) 



2ch 

where -D_a-i? -D-a-2 represent the parabolic cylinder function (or Weber function) and z 
and A are given as 



Tce-^H, (7) 
>? 
2c 



ih 2 

A = (8) 
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The final value of b(t) at t — is 

h^2- ih2 / 4c e~ nh2/8c 



6(0) = 



x 



2^ 

f 1 | y/3e 3 / 4 "(l+^ 2 /2c) ) 

lr(i + i/i 2 /4c) /ir(3/2 + ^ 2 /4c) y {) 



The probability to find the system in the ground state at t — is, when /i 2 /c ^> 1, 

Qr 2 

P QA (0) = |a(0)| 2 = l-|6(0)| 2 ~l- Mi . (10) 

Thus the probability -Pqa(^) does not approach 1 for finite c. 

We next present the solution for T(t) = c/t with i changing from to oo under the initial 
condition a = b = l/y/2 (see Appendix): 

b(t) = -^=e iht t ic F(l + ic,l + 2ic; -2iht), (11) 
v2 

where F is the confluent hypergeometric function. The asymptotic form of b(t) as t — > oo is 

v^(2fe)- fc r(2tc) 

— m — 

x { e -f**-W2 + ce ^+- c / 2 (2/it)- 1 } . (12) 
The probability to find the system in the target ground state behaves asymptotically as 

P QA (t) = \a(t)\ 2 (13) 

= l-|6(f)| 2 (14) 
sinh(7rc) 



1 - 



sinh(27rc) 



f nc ccos(2ht) cV c ) , . 

X { e + ^H + 4AV} (15) 
~ 1 - e- 2wc , (16) 

the last approximation being valid for c ^> 1 after t — > oo. The system does not reach the 
ground state as t — > oo as long as c is finite. Larger c gives more accurate approach to 
the ground state, which is reasonable because it takes longer time to reach a given value of 
T(= c/t) for larger c, implying slower annealing. 

The final example of the solvable model concerns the annealing schedule T(t) = c/y/i. 
The solution for b(t) is derived in Appendix under the initial condition a = b = l/y/2 as 

b(t) = -^f ht F Q - it, -2iht) + ^= e {3/4)m e iht (-2tht)^ 2 F (l - i 7 , | -2i/it) , (17) 
where 7 = c 2 /2h. The large-t behavior is found to be 
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b(t) ~ y/7ie 



-7TC 2 /4/l 



f 1 y/he (5/i) 
t (2htrU _ , , - + x 



^r(|-i 7 ) cT(-i 7 ) 



2 

f e -(i/4) 



1 ; \V2r(z 7 ) 2^r(i + ^ 7 )J. 

and the probability Pqa(oo) for c 2 //i ^> 1 is obtained as 



(18) 



P QA (oo) = l-|6(oo)| 2 ~l-^. (19) 

This equation indicates that the single-spin system does not reach the ground state under 
the present annealing schedule T(t) = c/y/i for which the numerical data in the previous 
section suggested an accurate approach. We therefore conclude that the asymptotic value 
of -Pqa(£) m the previous section may not be exactly equal to 1 for T(t) = 3/y/i although it 
is very close to 1. 

The annealing schedule T(t) = cj y/i has a feature which distinguishes this function from 
the other ones — ct and c/t. As we saw in the previous discussion, the final asymptotic 
value of -Pqa(^) is not 1 if the initial condition corresponds to the ground state for T — > oo, 
a = b = l/y/2. However, as shown in Appendix, by an appropriate choice of the initial 
condition, it is possible to drive the system to the ground state if T(t) = c/y/i. This is not 
possible for any initial conditions in the case of T(t) = —ct or c/t. 



V. SUMMARY AND DISCUSSIONS 

We have proposed the idea of quantum annealing (QA) in which quantum tunneling 
effects cause transitions between states in optimization problem, in contrast to the usual 
thermal transitions in simulated annealing (SA). The idea was tested in the transverse Ising 
model obeying the time- dependent Schrodinger equation. The transverse field term was 
controlled so that the system approaches the ground state. The numerical results on the 
probability to find the system in the ground state were compared with the corresponding 
probability derived from the numerical solution of the master equation representing the SA 
processes. We have found that QA shows convergence to the optimal (ground) state with 
larger probability than SA in all cases if the same annealing schedule is used. The system 
approaches the ground state rather accurately in QA for the annealing schedule r = c/y/i 
but not for faster decrease of the transverse field. 

We have also solved the single-spin model exactly for QA in the cases of T(t) = —ct, c/t 
and c/ y/i. The results showed that the ground state is not reached perfectly for all these 
annealing schedules. Therefore the asymptotic values of Pqa(£) in numerical calculations 
are probably not exactly 1 although they seem to be quite close to the optimal value 1. 

The rate of approach to the asymptotic value close to 1, 1 — Pqa(£), was found to be 
proportional to \/t in Fig. ^| for the ferromagnetic model. On the other hand, the single- 



spin solution shows the existence of a term proportional to 1/y/t, see Eq. (18). Probably 
the coefficient of the l/vt-term is very small in the situation of Fig. |] and the next-order 
contribution dominates in the time region shown in Fig. H. 
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A simple argument using perturbation theory yields useful information about the asymp- 
totic form of the probability function if we assume that the system follows quasi-static states 
during dynamical processes. The probability to find the system in the ground state is ex- 
pressed using the perturbation in terms of r(<^C 1) as 

p QA (r) ~i-r 2 £ _ (0) ! J0V (20) 



where E^ is the energy of the ith state of the non-perturbed (classical) system and E^ is 
the ground-state energy. If we set T — c/ we have 

Thus the approach to the asymptotic value is proportional to 1/t as long as the system stays 
in quasi-static states. The corresponding probability for SA is 

Psa(T) ~ ~ 1 - £ e-^l\ (22) 

which shows absence of universal (1/t-like) dependence on time. 

The present method of QA bears some similarity to the approach by the generalized 
transition probability in which the dynamics is described by the master equation but the 
transition probability has power-law dependence on the temperature in contrast to the usual 
exponential form of the Boltzmann factor This power-law dependence on the temper- 
ature allows the system to search for a wider region in the phase space because of larger 
probabilities of transition to higher-energy states at a given T(t), which may be the reason 
of faster convergence to the optimal states HH]. The transverse field term T in our QA rep- 
resents the rate of transition between states which is larger than the transition rate in SA 
(see (j5j)) at a given small value of the control parameter T(t) = T(t). This larger transition 
probability may lead to a more active search in wider regions of the phase space, leading to 
better convergence similarly to the case of the generalized transition probability. 

We have solved the Schrodinger equation and the master equation directly by numerical 
methods for the purpose of comparison of QA and SA. This method faces difficulties for 
larger iV because the number of states increases exponentially as 2^. The classical SA 
solves this problem by exploiting stochastic processes, Monte Carlo simulations, which have 
the computational complexity growing as a power of N. The corresponding reduction of the 
computational complexity is lacking in QA, and it is an important future problem in practical 
implementation of the idea of QA. Another future problem is to devise implementations of 
QA in other optimization problems such as the travelling salesman problem or the graph 
bipartitioning for which there seems to be no direct analogue of the transverse field to cause 
quantum transitions. 
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APPENDIX: SINGLE-SPIN PROBLEM 



In this Appendix we explain some technical aspects to derive the exact solution of the 
time-dependent Schrodinger equation for the transverse Ising model with a single spin. The 
three cases of T(t) = —ct, c/t and c/^/t will be discussed. 



1. Case of r(i) = -ct (Landau-Zener model) |g-[l3[ 

Let us express the solution of the Schrodinger equation at time t by the parameters 
a = (+\4>(t)) and b = ( — \ip(t)). The Schrodinger equation (|3|) with H = —ho z — Ta x is 
expressed as a set of first order differential equations for a and b. It is convenient to change 
the variables as 

a=^=(a + b), &=_( a -&), (Al) 



by which the Schrodinger equation is now 

-ic + h 2 + c 2 t 2 )b{t) = 0. (A2) 



dt 2 



By using the notation 



z = y/2ce~^H, (A3) 
ih 2 



we find 



^ + (A + i-i^( t ) = 0. (A5, 

The initial state is specified as a = b = l/y/2 or b = as t — > — oo. The solution of 
satisfying this condition is the parabolic cylinder function D_a-i(— zz) JT4j. Thus, we 
obtain the solution as 

a(t) = i f-rf6(t) - , (A6) 



/i y dt 
b(t) = C.D.^tz), (A7) 

where G\ is a constant. To fix C\, we use the condition 

5(-oo = 1 =1. (A8) 

hy 2c 

Then we have 

d = ^e~^ 2 / 8c . (A9) 



The wave function of this system is given in Eq. 



10 



2. Case of T(t) = c/t 

We next consider the case of T(t) = c/t. By eliminating a from the Schrddinger equation, 
we obtain 

d 2 b(t) _ 1 dT(t)db(t) 
dt 2 T(t) dt dt 

^^-wv^Y 1 ^*- (Aio) 

Substituting T(t) = c/t, we have 

d 2 b(t) Idbit) L ih c 2 \ I/x /A s 



The solutions of this equation are expressed by the confluent P function |L4 

f oo 

P { ih 1 ic t 
[ — ih — ic 

( oo 

= c M rP ' o i+ic -2zM !> , (A12) 
^ 1 ic —He 

the right-hand side of which has two independent expressions in terms of the confluent 
hypergeometric function 

j{t) = e iht t ic F{l + ic,l + 2ic; -2iht), (A13) 
g {t) = e iht t ic (-2iht)- 2ic 

xF(l-ic,l-2ic;-2iht). (A14) 

The general solution is b(t) = C\f{t) + C>2g(t). Using the initial condition 

6(0) = Ci + C 2 = -^, (A15) 

a(0) =C 1 -C 2 = ± (A16) 



we find 



b(t) = -^-=e iht t ic F{l + ic, 1 + 2ic; -2%ht). (A17) 
v2 



The asymptotic forms of bit) and |&(t)| 2 are then given as 

^/2{2h)- ic T{2ic) 



b(t) 



r(»c) 

x { e -iW-W2 + } , (A18) 
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3. Case of T(t) = c/<ft 



The final solvable model has T(t) = c/y/i. The Schrodinger equation ( |A10| ) is then 
expressed as 



d 2 b(t) 1 db(t) 
~dT 2 2t dt 



/ 9 2c 2 + z'/T 

+ [ h+ ^r- 



b(t) = 0. 



(A20) 



The solution is the confluent P function [14| 

oo 



Pi 



ih | — 27 t 
—ih 27 I 







jhtp j - Q 1 - Q 



(A21) 



1 z 7 5 

where 7 = c 2 /2h. The two independent solutions are thus [jl4 



1 



/(f) = e ^F(--z7,-;-2zM), 
g (t) = e lh \-2iht) l / 2 F(l - z 7 , -; -2i/it). 



(A22) 
(A23) 



The general solution of ( A20 ) is therefore the linear combination of the above two functions 

b(t) = C 1 f(t)+C 2 g(t). (A24) 
The constants C\ and C2 are fixed by the requirement 

(A25) 
(A26) 



6(0) = d 



1 

1? 



a (0) = 20. e W^c 3 = 4=- 



Substituting Ci and C 2 into Eq. ( |A24j) , we find 

c 
7^ 



6(f) = -j=e iW F (~ - z 7 , i; -2iw) + ^e^e M (-2iht)^F (l - z 7 , | -2iW ) . (A27) 



v/2 

The asymptotic form is 



-7TC 2 /4/l 



e^ ht (2/it)- 



-17 



+ 



v^r(i - 27) cr(-z 7 ) 



+e iht (2ht) 



-1/2+17 J !L 



-(l/4)7ri 



+ 



V5r(i T ) 2v^r(i + z 7 ) 



(A28) 
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The probability |6(oo)| 2 that the system remains in the excited state can be calculated as 
the asymptotic form of ( |A28| ) with the condition c 2 /h 3> 1 



|6(oo) 



ire 



"77T 



+ 



7 



-1/2 (5/4)™ 



"77T 



r(i-i 7 ) r(-i T ) 

«7 



,1/2 f'l 



"77T 



(-17) 



*7 



+ 7 - 1 /2 e (5/4) 7 ri^_^y 7 +l/2 



/> 2 



47 



2567 64c 4 ' 



(A29) 
(A30) 
(A31) 



4. Dependence of the final value on the initial condition 



We show that we can choose the initial condition so that the final state is the ground 
state when Y = c/y/t. This is not possible for T = —ct or c/t. From Eq. (|A24 ), the 
asymptotic form of the solution as t — > oo is 



b(t) ~ d 



-7rc 2 /Ah—iht 



(2ht) 



-ic 2 /2h 



T(l/2 - ic 2 /2h) 
i^he" 7rc2 / 4h - iht (2ht)~ ic2 / 2h 



c 2 T{-ic 2 /2h) 



(A32) 



The coefficients C\ and C 2 are fixed under the conditions 6(oo) = and |a(0)| 2 + |6(0)|" 



as 



C 2 



1 + 



sinh(iTC 2 /h) 
2sinh 2 (vrc 2 /2/i) 
ic 2 T(-ic 2 /2h) 
hT{l/2-ic 2 /2h) 1- 



-1/2 



(A33) 
(A34) 



This solution is not the ground state of the Hamiltonian 7i(0). 

The reason why one cannot obtain such a solution for the other schedules (T = c/t, —ct) 
is the following: The general solution for T = c/t also has two coefficients, and the initial 
state is represented as the linear combination of two terms whose phases are indefinite: 



a(0) = Ctt 
6(0) = C x t 



- c 2 r 
+ c 2 r 



(A35) 
(A36) 



The lowest-energy state at t = corresponds to a(0) = 6(0) = l/y/2 (times an arbitrary 
phase factor), which is realized by choosing C 2 = in Eqs. ( |A35| ) and ( |A36| ). The indef- 
initeness of t lc as t — > is irrelevant because this is only the overall phase factor. Such a 
situation does not happen for other values of a(0) and 6(0), leading to a serious difficulty 
to determine the wave function at t — 0. Thus we cannot choose an initial condition other 
than a(0) = b(0) = l/y/2. A similar fact exists in the case of Y — —ct. 
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FIGURES 



FIG. 1. Time dependence of the overlaps P S a(*), Pqa (t) , Psa( t (*)) and P^ A (T(t)) of the 
ferromagnetic model with T(t) = T(t) = 3/log(i + 1). 

FIG. 2. Time dependence of the overlaps of the ferromagnetic model with T(t) = T(t) = 3/\/i. 

FIG. 3. Time dependence of 1 — Pqa(£) of the ferromagnetic model with T(t) = 3/\/t. The 
dotted line represents to guide the eye. 

FIG. 4. Time dependence of the overlaps of the ferromagnetic model with T(t) = T(t) = 3/t. 

FIG. 5. The frustrated model where the solid lines denote ferromagnetic interactions and the 
broken line is for an antiferromagnetic interaction. 

FIG. 6. Correlation functions of spins 3 and 6 in Fig. [| for the classical and quantum cases. 
In the classical model (full line) the correlation is shown as a function of temperature while the 
quantum case (dotted line) is regarded as a function of the transverse field. 

FIG. 7. Time dependence of the overlaps of the frustrated model under T(t) = T(t) = 3/\/t. 
Here r c and T c are the points where the correlation functions vanish in Fig. ||. 

FIG. 8. Time dependence of the overlaps for the SK model with T(t) = T(t) = 3/y/i. 
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